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I.  Introduction 

The  constituents  of  the  upper  atmosphere  undergo  dissociation  and 
ionization  under  the  influence  of  the  intense  solar  radiation  available  at 
high  altitudes.  These  processes  are  of  interest  "because  they  are  respon- 
sible for  the  layers  of  free  electrons  in  the  upper  atmosphere  which  are 
so  important  in  radio  wave  propagation.  They  are  also  of  interest  because 
they  are  responsible  for  the  conversion  of  radiant  energy  into  heat. 

The  processes  in  the  upper  atmosphere  are  so  little  understood  that 
it  appears  impossible  to  treat  the  dissociation  phenomena  with  complete  gener- 
ality. However,  it  seems  possible  to  set  up  models  for  some  of  the  processes 
which  are  useful  in  gaining  insight  into  the  nature  of  the  phenomena.  A  pro- 
cess which  appears  amenable  to  theoretical  treatment  is  the  dissociation  of 
oxygen  in  the  upper  atmosphere.  This  process  is  comparatively  simple,  and 
one  can  set  up  reasonable  models  which  have  internal  consistency  and  can  be 
used  as  a  basis  for  making  numerical  calculations. 

In  a  previous  paper  [l]  a  self-consistent  model  for  the  dissociation 
of  oxygen  was  set  up.  The  object  of  the  present  paper  is  to  set  up  another 
model  and  to  present  numerical  calculations  based  on  it.  Inasmuch  as  the 
model  of  the  present  paper  has  many  features  in  common  with  that  of  refer- 
ence 1,  we  shall  review  the  assumptions  upon  which  the  model  of  reference  1 
is  based  and  indicate  the  modification  which  is  made  in  the  present  paper. 
One  of  the  more  controversial  assumptions  will  be  dropped'  and  another  sub- 
stituted. 

1.  The  principal  dissociation  process  is 

0_(x3  ]T~  )  +  hV->  0  (3p)  +  oO"D)  (1) 

d  g      ^U 

where  to  is  the  energy  of  the  quantum  of  radiation  absorbed  from  the  sunlight. 

2.  Tbe  principal  recombination  processes  are  two-body  recombination 
processes  in  which  two  oxygen  atoms  recombine   to  form  an  oxygen  molecule,   the 
excess  energy  being  in  the  form  of  radiation.     A  typical  process  of  this  sort 
is 

o(3p)  +  o(h)~*  QJtJ-T*)  +  to,  -»  o0(x3  T"  )  ♦  to       (2) 

a'      2   * — •      1     2    ' — 
where  n     denotes  the  recombination  cross  section. 
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3.  Dissociative  equilibrium  exists,  so  the  rate  of  dissociation  equals 
the  rate  of  recombination  at  any  altitude. 

U.  The  rate  of  energy  absorbed  in  process  (l)  equals  the  rate  of  emission 
of  energy  in  process  (2).  This  is  the  assumption  which  will  be  changed. 

5.  The  time  required  to  establish  dissociative  equilibrium  is  less  than 
that  required  to  establish  diffusive  equilibrium  between  the  oxygen  atoms  and 
molecules. 

6.  There  is  diffusive  equilibrium  between  the  nitrogen  and  oxygen  of 
the  atmosphere. 

There  are,  in  addition,  various  assumptions  of  a  less  fundamental 
and  a  less  controversial  nature.  They  were  used  in  reference  1  and  will  be 
brought  up  again  in  the  present  paper. 

With  the  above  assumptions  it  was  shown  to  be  possible  to  set  up 
a  system  of  three  equations  from  which  one  can  solve  for  the  three  fundamental 
unknowns:   the  temperature  T,  the  number  density  of  oxygen  atoms  n  ,  and  the 
number  density  of  oxygen  molecules  n_,  all  as  functions  of  altitude  x,  provided 
certain  boundary  conditions  were  given.  These  boundary  conditions  were  taken 
as  the  temperature  and  temperature  gradient  at  some  altitude  x  within  the 
region  where  dissociation  takes  place.  The  most  striking  feature  of  the  re- 
sulting calculation  was  that  the  dissociation  region  was  quite  large}  in  contrast 
to  the  results  of  previous  workers  who  found  that  the  region  was  quite  narrow 
(see  references  2-7 )• 

It  was  emphasized  in  reference  1  that  the  numerical  results  must 
not  be  taken  too  seriously  because  it  is  not  known  how  valid  the  assumptions 
1-6  are,  nor  how  accurately  the  cross  sections  (J  and  A'    represented  the  ex- 
perimental results.  In  particular^  had  to  be  assumed,  since  no  theoretical 
or  experimental  work  had  been  done  to  determine  this  cross  section.  The  prin- 
cipal result  of  reference  1  was  to  show  that  a  self-consistent  method  could  be 
set  up  with  the  assumptions  1-6,  and  hence  that  if  in  the  future  assumptions  1-6 
can  be  shown  valid  and  better  values  of  the  cross  sections  are  available,  the 
method  indicated  in  reference  1  can  be  used.  It  should  be  noted  that  the  assump- 
tions made  were  self-consistent;   it  was  not  necessary  to  assume  a  distribution 
n   which  is  one  of  the  quantities  to  be  calculated.   It  was  shown  that  the  work 
of  some  of  the  previous  workers  was  not  self-consistent  in  this  sense  . 

1.  Appendix  7,  Ref.  1 
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VI.  Other  Dissociation  Models 


Comparison  of  our  present  results  with  those  of  reference  1  show 

considerable  differences.   First  of  all,  in  reference  1  the  region  of  dis- 

12 
sociation  is  at  the  altitude  for  which  n_  =  10  .  Secondly,  the  dissocia- 
tion is  much  "broader,  of  the  order  15  Km.  Thirdly,  N„  is  much  larger  than 
the  value  obtained  using  the  ordinary  barometric  equation,  whereas  in  the 
present  case  it  is  smaller.   The  principal  reasons  for  these  differences 
are  that  we  have  dropped  the  energy  balance  requirement,  assumption  h,   and 
perhaps  more  important,  that  we  have  used  different  boundary  conditions. 
In  the  present  paper  we  have  required  explicitly  that  N  (x)— >  0  as  x— >oo. 
In  reference  1  we  could  not  apply  this  condition  explicitly  because  the 
theory  of  that  paper  was  not  valid  at  very  high  altitudes  since  the  energy- 
balance  condition  could  not  be  satisfied  there.  It  was  necessary  in  the 
theory  of  reference  1  to  assume  that  processes  occurred  such  that  n^(x) 
did  not  decrease  too  rapidly  at  high  altitudes,  but  the  discussion  of  such 
processes  were  outside  the  theory. 

The  results  of  the  present  model  are  more  satisfactory  than  those 
of  reference  1  in  that  N2(x)— >  0  for  x— >oo  as  required.   However,  this  does 
not  necessarily  mean  that  the  model  of  reference  1  is  incorrect  (although 
the  numerical  results  seem  untenable),  since  with  more  accurate  cross  sec- 
tions the  model  of  reference  1  might  give  results  closer  to  those  obtained 
in  the  present  work,  and,  specifically,  might  lead  to  smaller  values  of  N  . 

An  alternative  way  of  trying  to  modify  the  model  of  reference  1 
is  the .following: 

If  in  the  model  of  reference  1  one  had  taken  as  boundary  conditions 
N(x)— >  0  as  x— >oo  and  had  given  a  value  of  n0  at  some  altitude  x  ,  one  would 
have  found  that  the  temperatures  calculated  from  the  model  would  be  of  the  order 
of  a  thousand  degrees.  This  result  suggests  that  in  order  to  obtain  more  reason- 
able temperatures  one  should,  in  formulating  the  energy  balance  equations  take 
into  account  possible  ways  of  conserving  energy  other  than  by  the  radiation 
processes  (1)  and  (2);   in  particular,  one  should  include  heat  conduction. 
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The  present  calculations  indicate  that  the  dissociation  region 
occurs  where  the  atmospheric  density  is  so  great  that  it  is  conceivable 
that   three-body  collisions  are  important  as  a  recombination  process.     It 
is  possible  within  the  framework  of  the  present  paper     to  assume  three- 
body  recombinations  rather  than  two-body  recombinations. 
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ABSTRACT 

In  a  previous  paper  (ref.  l)  a  self-consistent 
model  was  set  up  for  the  dissociation  of  oxygen  in  the 
upper  atmosphere  under  the  assumption  of  number  balance 
and  radiative  energy  balance  and  a  form  of  the  "barometric 
equation.  In  the  present  paper  a  similar  model  is  used, 
but  the  requirement  of  radiative  energy  balance  is  dropped 
and  a  temperature  distribution  is  assumed.  The  resulting 
calculations  indicate  that  for  this  model  the  dissociation 
region  is  much  narrower  than  in  the  previous  model  and 
that  it  occurs  at  lower  altitudes. 
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I.  Introduction 

The  constituents  of  the  upper  atmosphere  undergo  dissociation  and 
ionization  under  the  influence  of  the  intense  solar  radiation  available  at 
high  altitudes.  These  processes  are  of  interest  because  they  are  respon- 
sible for  the  layers  of  free  electrons  in  the  upper  atmosphere  which  are 
so  important  in  radio  wave  propagation.  They  are  also  of  interest  because 
they  are  responsible  for  the  conversion  of  radiant  energy  into  heat. 

The  processes  in  the  upper  atmosphere  are  so  little  understood  that 
it  appears  impossible  to  treat  the  dissociation  phenomena  with  complete  gener- 
ality. However,  it  seems  possible  to  set  up  models  for  some  of  the  processes 
which  are  useful  in  gaining  insight  into  the  nature  of  the  phenomena.  A  pro- 
cess which  appears  amenable  to  theoretical  treatment  is  the  dissociation  of 
oxygen  in  the  upper  atmosphere.   This  process  is  comparatively  simple,  and 
one  can  eet  up  reasonable  models  which  have  internal  consistency  and  can  be 
used  as  a  basis  for  making  numerical  calculations. 

In  a  previous  paper[\ja  self-consistent  model  for  the  dissociation 
of  oxygen  was  set  up.   The  object  of  the  present  paper  is  to  set  up  another 
model  and  to  present  numerical  calculations  based  on  it.  Inasmuch  as  the 
model  of  the  present  paper  has  many  features  in  common  with  that  of  refer- 
ence 1,  we  shall  review  the  assumptions  upon  which  the  model  of  reference  1 
is  based  and  indicate  the  modification  which  is  made  in  the  present  paper. 
One  of  the  more  controversial  assumptions  will  be  dropped  and  another  sub- 
stituted. 

1.  The  principal  dissociation  process  is 

0-(x3  T~)  *  hl>-»  0  (3p)  ♦  O^D)  (1) 

2      g     *V 

where  hp  is  the  energy  of  the  quantum  of  radiation  absorbed  from  the  sunlight. 

2.  The  principal  recombination  processes  are  two-body     recombination 
processes  in  which  two  oxygen  atoms   recombine  to   form  an  oxygen  molecule,    the 
excess   energy  being  in  the  form  of  radiation.     A  typical  process  of  this   sort 
is 

o(3p)  +  o(h)  — >    o_(bi  5~+)  +  hv.  — »  o.Cx3  ]TI)  ♦  fc»       (2) 
fi'v  s  s 

where  Si    denotes  the  recombination  cross  section. 
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3.  Dissociative  equilibrium  exists,  so  the  rate  of  dissociation  equals 
the  rate  of  recombination  at  any  altitude. 

h.  The  rate  of  energy  absorbed  in  process  (l)  equals  the  rate  of  emission 
of  energy  in  process  (2).  This  is  the  assumption  which  will  be  changed. 

5.  The  time  required  to  establish  dissociative  equilibrium  is  less  than 
that  required  to  establish  diffusive  equilibrium  between  the  oxygen  atoms  and 
molecules. 

6.  There  is  diffusive  equilibrium  between  the  nitrogen  and  oxygen  of 
the  atmosphere. 

There  are,  in  addition,  various  assumptions  of  a  less  fundamental 
and  a  less  controversial  nature.  They  were  used  in  reference  1  and  will  be 
brought  up  again  in  the  present  paper. 

With  the  above  assumptions  it  was  shown  to  be  possible  to  set  up 
a  system  of  three  equations  from  which  one  can  solve  for  the  three  fundamental 
unknowns:  the  temperature  T,  the  number  density  of  oxygen  atoms  11  ,  and  the 
number  density  of  oxygen  molecules  n?,  all  as  functions  of  altitude  x,  provided 
certain  boundary  conditions  were  given.  These  boundary  conditions  were  taken 
as  the  temperature  and  temperature  gradient  at  some  altitude  x  within  the 
region  where  dissociation  takes  place.  The  most  striking  feature  of  the  re- 
sulting calculation  was  that  the  dissociation  region  was  quite  large  in  contrast 
to  the  results  of  previous  workers  who  found  that  the  region  was  quite  narrow 
(see  references  2-7). 

It  was  emphasized  in  reference  1  that  the  numerical  results  must 
not  be  taken  too  seriously  because  it  is  not  known  how  valid  the  assumptions 
1-6  are,  nor  how  accurately  the  cross  sections  G      and  fi'   represented  the  ex- 
perimental results.   In  particular  s'   had  to  be  assumed,  since  no  theoretical 
or  experimental  work  had  been  done  to  determine  this  cross  section.  The  prin- 
cipal result  of  reference  1  was  to  show  that  a  self-consistent  method  could  be 
set  up  with  the  assumptions  1-6,  and  hence  that  if  the  future  assumptions  1-6 
can  be  shown  valid  and  better  values  of  the  cross  sections  are  available,  the 
method  indicated  in  reference  1  can  be  used.  It  should  be  noted  that  the  assump- 
tions made  were  self-consistent:   it  wa6  not  necessary  to  assume  a  distribution 
np,  which  is  one  of  the  quantities  to  be  calculated.  It  was  shown  that  the  work 
of  some  of  the  previous  workers  was  not  self-consistent  in  this  sense  . 

1.  Appendix  7,  Hef.  1 
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As  indicated  above,   perhaps  the  most  controversial  assumption  is  assump- 
tion k,   that  of  energy  conservation  through  a  radiative  steady  state.  This  assump- 
tion is  questionable  "because  energy  processes  other  than  (l)  or  (2)  may  occur, e.g. 
heat  conduction.  In  this  present  paper  we  shall  try  to  eliminate  assumption  U. 
It  will  now  no  longer  he  possible  to  calculate  all  three  quantities  T,  n?,  n  as 
functions  of  altitude,  as  was  done  in  reference  1,  for  we  shall  have  too  few  equa- 
tions.  Instead  we  shall  take  the  temperature  T  as  given,  and  calculate  n  and  n 
as  functions  of  the  altitude  under  appropriate  boundary  conditions. 

Dropping  assumption  U,  while  having  the  disadvantage  of  compelling  us  to 
assume  a  value  for  T,  has  the  advantage  of  making  the  resulting  theory  applicable 
over  a  larger  range  of  altitudes.   This  is  because  assumption  H  can  be  true  only 
in  regions  where  processes  (l)  and  (2)  are  principal  processes  leeding  to  con- 
version of  radiation  to  heat,  whereas  the  treatment  of  the  present  paper  will  be 

true  at  any  altitude,  provided  the  atmospheric  density  is  sufficiently  low  to 

2 
permit  one  to  assume  two-body  recombinations.   As  a  result,  we  shall  be  able 

to  impose  boundary  conditions  on  our  equations  at  extremely  high  altitudes. 
It  is  to  be  emphasized  that  in  this  paper,  as  in  reference  1,  the 
numerical  calculations  should  not  be  taken  too  seriously.  The  numerical  cal- 
culations are  undertaken  essentially  to  shov  the  nature  of  the  solution  of  the 
problem. 

II.  The  Basic  Equations 

We  shall  now  write  the  basic  equations  in  some  generality  and  later 
include  our  assumptions.  The  first  of  these  equations  is  the  number  balance 
equation  in  which  the  number  of  dissociations  per  unit  volume  per  second  enuals 
the  number  of  recombinations  per  second.  Let  us  consider  the  dissociation 
processes.   It  is  generally  agreed  thet  the  principal  dissociation  process  is 
process  (l)  in  which  the  oxygen  molecule  is  dissociated  from  the  normal  elec- 
tronic bound  state  to  the  dissociated  state  in  the  Runge-Schumann  eontinuum  on 
the  absorption  of  radiation  with  a  wave  length  smaller  than  If^Ok.     The  cross 
section  for  this  process  has  been  studied  experimentally  in  references  8  and  9 

2.  See  concluding  paragraph  of  present  paper. 
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and  theoretically  in  reference  10.     Denoting  this  cross  section  by  (T    as  above, 
ve  have  the  number  of  dissociations  per  second  at  any  altitude  z  given  by 


Hw/  ff^^W* 


where  P   (i)  is  the  radiation  density  per  unit  frequency  available  at  the  alti- 

tude  x,  n?(i)  is  the  number  density  of  oxygen  molecules  at  the  altitude  x,  and 
v  is  the  threshold  frequency  corresponding  to  the  wave  length  If^Ok.     How, 

since  there  is  attenuation  of  light  by  the  oxygen  molecules  above  the  altitude 
x,  we  have 

/>„<x>  -/>,(»)•  V     2 


00 

Hg(x)  =  f   n2(x)  dx 


(3) 


where  f   ( oo  )  represents  the  radiation  density  available  at  the  "top"  of  the 
atmosphere.  Then  one  can  write 

no.  dissociations/sec  °  n«(x)  l(Hp)  (*0 

The  number  of  recombinations  per  unit  volume  per  second  can  be  written  as 

n  (x)  B'(T)  where  n  (x)  is  the  number  density  of  oxygen  atoms  at  altitude  x 
P  P 

and  B'(T)  it  the  recombination  coefficient,  which  depends  on  the  specific  re- 
combination assumed,  and  T  is  the  temperature  at  the  altitude  x.  The  recombin- 
ation coefficient  depends  upon  the  temperature  because  it  has  been  averaged  over 
the  relative  velocities  of  the  atoms  using  the  Soltsmann  distribution  character- 
ised by  a  temperature  I.  Hence  the  equation  of  number  balance  is 

2 


UjjU)  l(H2(x))-  sjB'(T).  (6) 
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The  barometric  equation  it 

where  P  is  the  pressure  at  altitude  z,  g  it  the  gravitational  constant,  N  is 
the  mass  of  the  ox7gen  atom,  n^  is  the  number  density  of  nitrogen  molecules, 
and  Mg  is  the  mass  of  the  nitrogen  molecule, 

Bat   P  -  (nj  +  lig  +  Bj)  H  .  (7) 

Hence  the  barometric  equation  is 

h  [*"*  +  *2  *  V  M ]  "  -  <*"*  *  «*2  *  Vh)  *•       (8> 

Iquations  (6)  and  (8)  are  the  basic  equations  of  our  theory.  It  should  be  noted 
that  these  equations  are  quite  general  and  are  largely  independent  of  the  model. 
In  these  equations  ve  shall  assume  T  and  a.  as  given  functions  of  x.  The  un- 
known functions  which  we  wish  to  so  Ire  for  are  n  and  n_  , 

III.  Detailed  Model 

Ve  shall  work  with  a  model  in  which  the  radiation  from  the  sun  is 
equivalent  to  radiation  from  a  black  body  of  temperature  6000°  K.  Together 
with  certain  geometrical  considerations,  (see  ref.  l)  this  assumption  speci- 
fies />p(oo).  A  more  controversial  assumption  is  that  which  determines  a. 

as  a  function  of  the  altitude  z.  Most  of  the  previous  workers  assume  that  the 

ma- 
ntle of  nitrogen  to  oxygen  7*.    has  a  constant  value,  namely  k.  the 

n2*2B* 
value  of  the  ratio  at  the  surface  of  the  earth.  This  assumption  requires 
that  winds  of  sufficient  frequency  and  strength  occur  so  as  to  disturb  dif- 
fusive equilibrium.  If  such  winds  occurred  they  would  most  likely  also  dis- 
turb the  distribution  of  oxygen  atoms  and  molecules  so  that  the  problem  of 
finding  the  dissociation  ratio  would  be  meaningless.  Ve  therefore  take  the 
alternative  point  of  view,  namely  that  diffusive  equilibrium  predominates  be- 
tween oxygen  and  nitrogen  (assumption  6).  This  assumption  means  that  oxygen 
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and  nitrogen  are  distributed  independently  of  each  other.  Vith  this  assumption 
we  hare 

fc  G*  4  -  -  %  *  •  <«> 

and  equation  (8)  becomes 

E  [(V  n2)  "]  "  "  [*V  ^  e  (8a) 

and  all  reference  to  the  distribution  of  nitrogen  in  the  barometric  equation  (8a) 
disappears. 

The  most  controversial  of  the  present  assumptions  is  that  which  deals  with 
the  recombination  process.  Some  of  the  previous  workers  (ref .  3»^»7)  assume  that 
recombination  of  the  two  oxygen  atoms  is  given  to  a  third  body  when  the  three 
bodies  collide.  Schematically 

0  +  0  +  M— >  0  ♦  M*  (10) 

where  M  is  the  third  body  and  the  asterisk  represents  the  third  body  in  a  vtate 
which  is  excited  either  internally  or  in  the  sense  that  it  has  greater  transla- 
te ional  energy. 

Others  (ref.  5,6)  prefer  to  think  of  the  recombination  as  being  a  two- 
body  radiative  process  in  which  the  excess  energy  is  emitted  in  the  form  of  radi- 
ation. 

0  +  0  — »02  +  h».  (11) 

For  simplicity  and  in  order  to  minimize  the  difference  between  the 
model  of  our  present  paper  and  that  of  reference  1  we  assume  two-body  recombina- 
tions in  the  present  calculation. 

If  one  assumes  a  two-body  process,  it  can  be  shown  (ref.  l)  that  this 
process  goes  so  slowly  that  virtually  all  the  oxygen  atoms  in  the  state  D,  which 
is  the  result  of  the  dissociation  (equation  (l)~)  ,  will  have  had  time  to  drop  to 
the  -*P  state  by  emission  of  the  forbidden  red  line 

O^D)  -»0(3P)  ♦  hPred  (12) 
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■o  that  Tirtually  all  the  atoms  are  in  the  ^P  state;  hence  the  notation  n     for 

the  atomic  number  density.     Then  the  most  probable   recombination  process,   assum- 
ing two-body  recombinations,   is  one  typified  by  the  equation 

0(3P)  ♦  0(3P)  (St   )— *0  (bl£  +)  +  hv  ,  (13) 


g 


and  the  rate  of  recombination  is 

SO 


2 

n 
P 


U T  ^  'I<T)  * ' 


where  t  is  the  relatiTe  Telocity  of  the  two  atoms,  A'     is  the  recombination 
cross  section  of  process  (13 )  and  f_(r)  is  the  Boltmann  distribution  with 
respect  to  relative  Telocity  r.  Unfortunately  the  cross  section  B'   for  pro- 
cess (13)  is  not  known,  but  one  can  estimate  it  by  assuming  that  within  a 
statistical  weight  factor  and  a  factor  arising  from  a  different  threshold  fre- 
quency it  is  much  the  same  as  the  known  cross  ssction  B     for  the  process  in- 
Terse  to  process  (2),  namely 

0(3P)  ♦  0(h)   — >  0  (X3  T~  ")  +  ho  .  (2a) 

g 
00 

Hence  we  can  eraluate  B'(T)  «y       rfiL  '■(▼J**-  of  equation  (6). 

To  replace  assumption  ik)  we  assume  instead  that  T  =  T     +<*(x-x   )  t 

0  00 

being  the  temperature  of  the  altitude  x  =  x     and  ex  being  the  gradient  of  the 

temperature.     In  point  of  fact  it  is   sufficient  to  take  T  =  T     =  const,     be- 

0        ' 

cause  it  turns  out  that  the  numerical  results  are  independent  of  ex  if  ex  i8 

sufficiently  small  (cx<  10  /Em,  say).  This  point  will  later  be  discussed  more 
fully.  Howerer,  we  shall  take  ex  arbitrary,  since  the  equations  are  not  essen- 
tially more  complicated  for  this  case. 

17.  The  Basic  Iquationa  and  the  Method  of  Solution 

The  basic  equations  for  the  calculations  are  the  number  balance  equa- 
tion 

l(H2)n2  -  B'(T)nJ  (6) 

and  the  barometric  equation  which  for  T  =  T     +«(x-x  )  has  the  following  form: 


d(n_  ♦  n_)  f 

dx  —  ■"     H^2+V+     k*(*p*2V 


(Sb) 


_  g  _ 


l(BL)  and  B'(T)  have  "been  calculated  and  are  given  graphically  as  functions  of 
their  arguments  in  figures  1  and  2. 

To  solve  the  equations  (6)  and  (9b)  it  is  convenient  to  introduce  the 
total  number  density 

f(x)  =  UgCx)  +  n^x).  (1*0 

From   (lU)  and  (6)  one  obtains 

/    /i(NP)       /un  ; 

Jimr  +7BtfT-  +  Uf 

n2  ~  V 2 /  (15) 


-y-^fi 


n_  = 


(16) 


2 


and  (Sb)  can  be  written  either  as 


^    =  -     £{(£♦«>£(*>♦£    n2(x)   >  (17) 


or 

^r1  "  r  |  (f^«)fU)-^*p(*)r.  us) 

It  is  convenient  to  consider  (17)  and  (15)  as  the  differential  equa- 
tion for  f(x)  when  n_(x)  is  small  and  (18)  and  (16)  as  the  differential  equation 

for  f(x)  when  n  (x)  is  small,  though  both  sets  of  equations  are,  of  course,  the 

-»Q0 

X 

At  high  altitudes  where  n_  is  very  small  compared  to  n  ,   it  is  seen 
that  f(x)  =  n     and  that  equation  (17)   is  the  ordinary  barometric  equation  for  n_. 
The  density  n_(x)  is  given  by  (15)  which,   since  f (x)  is  small  too   (but  much 


same.     Of  course     we  use  f  n„(x)  =  N„  in  either  case. 


larger  than  n2(x)),  becomes 


10 


-8a- 


1(0) 


10  - 


10 


I(N2) 


.67 


10 


10 


Figure  1 
Dissociation  Probability  as  a  Function  of 
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Figure  2 

Recombination  Cross-Section  as 
Function  of  Temperature. 
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n2(x)  "y  mr  f(x)  "V  ttc 


TTT7  np  •  (19) 


which  is  Just  equation  (6).  Thus  for  high  altitudes,  the  oxygen  atoms  "behave 
Just  like  an  ordinary  gas  without  dissociation  present,  and  n?  is  given  hy  equa- 
tion (6). 

For  high  altitudes  where  N_->  0  so  that  I(H_)  — >  1(0)  ■  const,  and  for 
the  case  that  T  is  a  constant  at  these  altitudes,  n  (x)  also  satisfies  the  ordin- 
ary barometric  equation  as  is  clear  from  (19). 

For  low  altitudes  where  IT  is  large  and  n  small,  f(x)  =  n„(x)  to  a 

very  good  approximation  and  equation  (IS)  becomes  the  ordinary  barometric  equa- 
tion for  n„(x).  That  is  if  n_(x)  were  prescribed  at  x  for  x  sufficiently  low. 

Z  2  0  0 

both  the  barometric   equation  with  and  without  dissociation  would  give  the  same 

result  for  n_(x)  at  low  altitudes.     n(»)  is  given  by  (17 )  which  beoomes 
(since     jiS   k     is  large) 

7  no        ■  V  -wry-  (20) 


n   *  /IS 


np 


which  Is  Just  equation  (6). 

Thus  for  low  altitudes  n?(x)  obeys  the  ordinary  barometric  equation 
whereas  for  high  altitudes  n (x)  obeys  it;   in  the  intermediate  region  neither 
does,  and  n  (i)  and  n?(x)  influence  each  other  greatly  with  the  consequence 

that  n?(x)  falls  much  more  rapidly  than  predicted  by  the  ordinary  barometric 
equation. 

Equations  (17)  and  (15)  or,  equivalently,  (lg)  and  (l6)  were  solved 
by  the  following  scheme: 

One  prescribes  n„(x  )  and  takes  for  H  (x  )  the  result  predicted  by 
the  ordinary  barometric  equation,  Just  to  have  some  starting  value  for  Hp(x). 

Using  equation  (6)  for  x  «  x  ,  n  (x  )  can  be  found,  and  from  (lh) 

one  obtains  f(x  ).  Then  using  (17 )  one  can  solve  for  f (x)  by  a  step-by-step 
0 

integration.  Tor  example,  for  x,  ■  .x  +A*  one  has 


Al* 


**xo 


where  ^ 
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is  computed  from  (17). 


x=x 


One  then  has 

I2(x1)  «  N2(xo)  -n2(xo)^  x, 

and  can  compute  np(x.)  from  (15),   since  f(x1)  and  I(B_)  at  x  =  x.    are  now  known. 

Further,   h^C^)  =  i(\)  -  ^^  and  dx  I  can  te  found  from  ^lg^* 

The  procedure  is  repeated  for  x?  =  x.    +ApX  and  one  integrates  to  as  high  an 
altitude  as  one  wishes. 

One  ultimately  gets  distributions  n  (x),  n_(x),  N_(x)  based  on  the 
choice  of  N„(x  ).     The  curre  for  IT_(x)  does  not  approach  0  as  x->oo.     Instead, 

2   0c 

one  gets  a  curve  which  approaches  a  constant  HL  , 


N. 


Ng  min 


Figure  3 

In  order  to  get  a  better  value  of  N_(x  )  one  takes  the  old  value  of 

N^(x  )  and  subtracts  H.  .  .  Using  the  new  value  of  H_(x  )  one  repeats  the  cal- 
2  o  2mm.  2  0 

culation.  Again  one  finds  a  curve  for  N_(x)  similar  to  the  one  above  except 

that  N„    is  much  smaller  (of  course,  one  wishes  it  to  be  0).  One  takes  still 
2  min 

a  third  value  of  N_(x  )  which  consists  of  the  second  value  of  H„(x  )  minus  the 
2   0  d      0 

second  value  of  N_    and  the  integration  is  repeated.  One  eventually  gets  a 
2  mm 

value  of  N0    which  is  so  small  and  occurs  at  such  a  large  value  of  x,  that  at 
2  min 

and  below  the  region  at  which  dissociation  becomes  important  the  results  are  no 
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longer  sensitive  to  N_  .   .  At  this  point  the  iteration  procedure  stops;   6 
iterations  were  required  in  the  calculations. 

The  final  result  is  the  following.  At  low  altitudes  the  curves  for 
n2(x)  with  and  without  dissociation  agree  if  the  curvee  are  fitted  at  some  point 

x'  (which  may  be  less  than  x  )  where  there  is  virtually  no  dissociation. 


Figure  U 


no  dissociation 

dissociation 


For  low  enough  altitudes  the  curves  for  H„(x)  with  and  without  dis- 
sociation will  also  agree  because  in  either  case  the  most  important  contribu- 
tion to  K2(x)  at  low  altitudes  (well  below  x  )  comes  from  the  value  of  xu(x) 
within  a  few  kilometers  of  x.   This  is  because  np(x)  in  both  cases  is  rela- 
tively very  small  for  x  at  any  distance  from  the  altitude  being  considered. 

The  important  feature  of  the  present  treatment  is  that  aside  from 
assigning  a  boundary  value  n~  at  x=x  ,  no  a  priori  assumption  is  made  on  the 
distribution  of  n„(x)  or  equivalently  on  the  function  N Ax).   These  distribu- 
tions are  obtained  from  the  basic  equations.   In  all  of  the  previous  work  on 
the  oxygen  dissociation  problem  (except  in  reference  l)  such  a  priori  assump- 
tions are  made,  as  pointed  out  in  Appendix  7  of  reference  1.   For  example, 
Penndorf  (ref .  7  )  appears  to  assume  that  one  may  take  U9(x)  (equivalently 
np(x))  as  though  it  were  given  by  the  ordinary  barometric  equation  when  no 
dissociation  is  present;   i.e.,  he  assumes  dissociation  has  only  a  minor 
effect  on  the  overall-distribution  of  n_.  It  will  be  seen  below  that  this  is 
not  true,  at  least  with  the  cross  sections  used. 
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V.  Results  of  Numerical  Calculations 

Three  cases  were  calculated  numerically.   In  all  three  cases  ru  at  x=x 

was  taken  as  10  /cm  .   (If  02  is  assumed  to  constitute  l/5  of  the  atmosphere, 

x  would  correspond  to  85  km.)  This  choice  is  not  special  though  it  may  seem  so 

since  n?(x),  nCx)  depend  only  on  the  difference  x-x  .  The  three  different 

cases  correspond  to  different  choices  of  the  temperature  function  T  =  T  +°<(x-x). 

In  the  first  case  T  was  taken  as  195  and  <*  was  taken  as  5  /l™.   In  the  second 

case  T  =  200°k,cx  =  0. 
o 

In  the  third  case  T  =  300°k,«  =  0.   Comparison  of  results  for  the 
first  and  second  cases  gave  essentially  the  same  results,  indicating  that 
small  temperature  gradients  do  not  affect  the  results  greatly.  This  conse- 
quency  is  due  to  the  smallness  of  the  region  in  which  dissociation  is  impor- 
tant;  if  the  temperature  has  not  too  large  a  gradient  in  this  small  region, 
the  temperature  is  essentially  constant. 

Since  the  first  and  second  cases  agree  so  closely,  we  have  given 
graphs  only  of  the  second  and  third  cases.   Figures  5  and  6  represent  the 
results  for  the  second  case,  and  figures  7  and  8  the  results  for  the  third  case. 

As  is  seen  in  "both  cases,  the  dissociation  takes  place  in  a  very  narrow 

lU 
region  of  the  order  of  a  kilometer  near  the  altitude  at  which  n0  =   10  ,  (hence 

the  choice  n_(x  )  =  10  ).     The  effect  of  the  higher  temperature  is  to  spread 

out  the  region  of  dissociation  and  raise  the  altitude  of  this  region.  N?  Is 

much  larger  for  the  higher  temperature,  the  reason  being  that  n_(x)  decreases 

more  slowly  with  the  higher  temperature. 

Also  indicated  on  the  graphs  are  the  results  when  the  ordinary  baro- 


lX-XQ) 

metric  equation  without  dissociation  is  used  for  ttp^iu  =  np0e  where 

n'  and  x'  are  chosen  so  that  they  coincide  with  the  value  obtained  for  n~(x)  at 
20      0  e. 

sufficiently  low  altitudes.   It  is  clear  that  in  the  region  of  dissociation  and 
above  that  region,  iu(x),  Np(x)  deviate  considerably  from  the  values  obtained 
from  the  ordinary  barometric  equations.  Hence  the  use  of  the  ordinary  baro- 
metric equation  by  some  workers  to  give  values  of  ruCx) ,  Np(x)  is  incorrect. 
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71.  Other  Dissociation  Models 


Comparison  of  oar  present  results  with  those  of  reference  1  show 

considerahle  differences.  First  of  all,  in  reference  1  the  region  of  dis- 

12 
sociation  is  at  the  altitude  for  which  n-  =  10   .  Secondly,  the  dissocia- 
tion is  much  broader,  of  the  order  15  Kin.  Thirdly,  Np  is  much  larger  than 
the  value  obtained  using  the  ordinary  barometric  equation,  whereas  in  the 
present  case  it  is  smaller.  The  principal  reasons  for  these  differences 
are  that  we  have  dropped  the  energy  balance  requirement,  assumption  h,   «T>d 
perhaps  more  important,  that  we  have  used  different  boundary  conditions. 
In  the  present  paper  we  have  required  explicitly  that  5 Ax)  —  >  0  as  x  —  >   00. 
In  reference  1  we  could  not  apply  this  condition  explicitly  because  the 
theory  of  that  paper  was  not  valid  at  very  high  altitudes  since  the  energy 
balance  condition  could  not  be  satisfied  there.   It  was  necessary  in  the 
theory  of  reference  1  to  assume  that  processes  occurred  such  that  n_(x) 
did  not  decrease  too  rpaidly  at  high  altitudes,  but  the  discussion  of  such 
processes  were  outside  the  theory. 

The  results  of  the  present  model  are  more  satisfactory  than  those 
of  reference  1  in  that  Np(x)— »0  for  x— ->oo  as  required.  However,  this  does 
not  necessarily  mean  that  the  model  of  reference  1  is  incorrect  (although 
the  numerical  results  seem  untenable),  since  with  more  accurate  cross  sec- 
tions the  model  of  reference  1  might  give  results  closer  to  those  obtained 
in  the  present  work,  and,  specific ally, might  lead  to  smaller  values  of  N„. 

An  alternative  way  of  trying  to  modify  the  model  of  reference  1 
is  the  following: 

If  in  the  model  of  reference  1  one  had  taken  as  boundary  conditions 
H(x) — »  0  as  x — >oo  and  had  given  a  value  of  n„  at  some  altitude  x  ,  one  would 
have  found  that  the  temperatures  calculated  from  the  model  would  be  of  the  order 
of  a  thousand  degrees.  This  result  suggests  that  in  order  to  obtain  more  reason- 
able temperatures  one  should,  in  formulating  the  energy  balance  equations,  take 
into  account  possible  ways  of  conserving  energy  other  than  by  the  radiation 
processes  (l)  and  (2);   in  particular,  one  should  include  heat  conduction. 
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The  present  calculations  indicate  that  the  dissociation  region 
occurs  where  the  atmospheric  density  is  so  great  that  it  is   conceivable 
that   three-body  collisions  are  important  as  a  recombination  process.      It 
is  possible  within  the  framework  of  the  present  paper,    to  assume  three- 
body  recombinations  rather  than  two-body  recombinations.      Such  a  model 
would  require  more  secondary  assumptions  such  as  the  efficiencies  of  the 
various  possible  third  bodies  and  their  distribution.     A  calculation, 
similar  to  the  present   one  but  using  the  assumption  of  three-body  recom- 
binations,   is  being  carried  out  and  will  be  reported  separately. 
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